Domain-wall dynamics in translationally non-invariant nanowires: theory and applications 
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We generalize domain-wall dynamics to the case of translationally non-invariant ferromagnetic nanowires. 
The obtained equations of motion make the description of the domain-wall propagation more realistic by ac- 
counting for the variations along the wire, such as disorder or change in the wire shape. We show that the 
effective equations of motion are very general and do not depend on the model details. As an example of their 
use, we consider an hourglass-shaped nanostrip in detail. A transverse domain wall is trapped in the middle and 
has two stable magnetization directions. We study the switching between the two directions by short current 
pulses. We obtain the exact time dependence of the current pulses required to switch the magnetization with the 
minimal Ohmic losses per switching. Furthermore, we find how the Ohmic losses per switching depend on the 
switching time for the optimal current pulse. As a result, we show that as a magnetic memory this nanodevice 
may be 10 5 times more energy efficient than the best modern devices. 

PACS numbers: 75.78.Fg, 75.60.Ch, 85.75.-d 
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During last two decades there has been a significant 
progress in describing magnetization and, in particular, 
domain- wall (DW) dynamics in magnetic nanostructures [1- 
3511 . The interest to these studies has been inspired not only 
by the fundamental physics questions but also by the potential 
applications for the spintronic memory and logic nanodevices 
1 36-38]. However, recently this progress has been staggered 
due to the inability to make perfect translationally invariant 
nanowires from one side and the difficulty to apply the theo- 
ries made for translationally invariant systems to successfully 
describe some of the phenomena in the experimental systems 
from the other side. There have been attempts to consider 
analytically and numerically 127143 ill the systems with rough 
surfaces or other disorder but the general theory for transla- 
tionally non-invariant magnetic systems is still lacking. 

In this Letter we generalize current and field induced DW 
dynamics to the case of translationally non-invariant ferro- 
magnetic nanowires. This generalization makes the descrip- 
tion of the DW motion more realistic by accounting for vari- 
ations along the wire, such as disorder or change in the wire 
shape. We show that the effective equations of DW motion 
are very general and do not depend on the details of the model 
Hamiltonian. These equations are the main theoretical result 
of this Letter. 

As an example of the application of this theory, we study 
current-induced magnetization switching in a thin hourglass- 
shaped nanostrip. We show that a transverse DW trapped in 
the middle of a curved inward nanostrip can serve as a mag- 
netic memory device, see Fig.Q] At zero current, the magne- 
tization in the transverse DW can have either of the two equi- 
librium directions in the plane of the nanostrip. We study the 
switching between these two magnetization directions medi- 
ated by short current pulses with the requirement that the DW 
returns to the initial position after switching. The main en- 
ergy loss during the switching is due to Ohmic heating of the 
wire. The minimum energy required per switch depends on 
the switching time. We obtain the exact time dependence of 
the current pulses required to switch the magnetization in the 



most efficient way (with the minimal Ohmic losses per switch- 
ing for a given switching time). Furthermore, we find how the 
Ohmic losses depend on the switching time for this optimal 
current pulse. 

The two equilibrium magnetization directions can serve as 
"zero" and "one" of a memory bit. We show that based on 
this prototype, it is possible to design a nonvolatile memory 
device with an extremely short writing time, which is only 
limited by the spin-wave frequency. The energy required per 
switch for such a nanodevice is much lower than that for the 
state-of-the-art memory devices. 

Equations of motion. We study the DW dynamics by em- 
ploying the Landau-Lifshitz-Gilbert (LLG) equation with cur- 
rent terms |T^, 28]. In general, a nontrivial solution of the 
static LLG equation, So x 5Ho / <5So = 0, has a continuous 
symmetry. It means that this solution can be parametrized by 
an even-dimensional vector £ such that So (z, £) is the solution 
for any £ in a continuous interval. Then the DW dynamics 
due to a correction to Hamiltonian T~Lq or an electric current 
can be described by the time-dependent parameter £(t). The 
equations for £(t) are called the effective equations of mo- 
tion. Below we sketch the general derivation of such effective 
equations. 

To consider the dynamics we assume that there is a pertur- 
bation h, such that the full LLG equation takes the form 
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where the time is measured in units of the gyromagnetic ra- 
tio 70 = g\e\/(2mc) and S = M/M with M being the 
saturation magnetization. We note that h is not only a cor- 
rection to Hamiltonian and may also contain other contribu- 
tions such as dissipation, adiabatic and nonadiabatic current 
terms. We look for a solution of this equation in the form 
S(z, t) = So {z, + s, where the dependence £(t) is weak 
and s is small and orthogonal to So at each point. 

For current-driven magnetization dynamics the current is 
assumed to be uniform in the nanowire. The correction h then 
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FIG. 1. (color online) Hourglass-shaped nanostrip as a prototype of 
a magnetic memory nanodevice. 



can be written as a sum of two terms: 



S x— - + h Q , 
dS 



, = aSoxSo-jdSo+/3jS o x0So. 

(2) 

The first term is the correction to the Hamiltonian, Hs, which 
turns the zero modes into soft modes as well as couples the 
magnetization to an external magnetic field. The second term, 
h Q , is due to dissipation and current terms. Here a is the 
Gilbert damping constant, j is an electric current in the units 
of velocity, and (3 is the non-adiabatic spin torque constant. 

Taking the scalar product of Eq. © with So (z) x dt.So{z) 
and integrating over the space, one can find the following 
equation for the collective coordinates £j : 



2& = '"'k I- - e ij fdzS (z) x %S (z) • h a (z), 



(3) 



where E(£) = Hs[Sq(z, £)] is the energy of the domain wall 
ll39ll as a function of the soft modes 

For thin nanowires, a DW is a rigid spin texture. Its slow 
dynamics can be described in terms of only two collective co- 
ordinates corresponding to zero modes of motion. These zero 
modes are the DW position zq and its conjugate variable - the 
tilt angle <f> for the transverse DW. For the vortex DW, cf> can 
serve as the magnetization angle defining the transverse posi- 
tion of the vortex in the wire B40I1 . Using the definition of h Q , 
Eq. ©, and the fact that So = — zq8 z So + (pd^So we obtain 



h Q = a^Soxd^So-j^So + ^j-ai^SoXc^So-SoxH 

(4) 

Here we included uniform magnetic field H in h a by adding 
the term —So x H in ©. 

Up to the leading order in small dissipation (a and j3) the 
equations of motion become 



1 dE 

1 dE aa zz dE 



2 dz 2 dcp 



(5) 



+ H + (a - /3)a zzJ . (6) 



Here for simplicity we have taken H to be along the wire (in 
the z direction) and a zz = 5 J dz(d z So) 2 . These equations 



are rather universal and can also be applied to describe the dy- 
namics of vortex domain walls [41] in terms of two collective 
coordinates associated with the DW degrees of freedom. The 
1/2 in the first term on the right-hand side of Eqs. © and © is 
a consequence of the Poisson bracket of the conjugated vari- 
ables zq and <fi. The most general derivation of these equations 
is based on Poisson brackets and energy dissipation and will 
be presented elsewhere j42Tl . 

Equations © and © do not depend on details of the micro- 
scopic model. The only required input is the energy of a static 
DW as a function of two parameters z and <fi. This function 
can be either calculated by means of an analytical approxi- 
mate model, micromagnetic simulations, or can be measured 
experimentally for a give n wire by a method analogous to the 
one described in Ref.l39l 

The energy of a static DW, E(zQ,(f>) = H[Sq(z; zn, </>)], 
where So is a solution of a static LLG, in general depends on 
both zo and <p>. The main dependence of E on the angle cf> 
comes from the anisotropy in the transverse plane, E(<ft) = 
—KCOs(2(j)) B43I1 . The dependence of E on the DW position 
zq may come from different sources such as, z-dependence of 
the wire shape, nonuniform concentration of impurities, wire 
surface roughness, nanofabricated notches, etc. Equations © 
and © can be used to study DW propagation in disordered 
wires as well as DW depinning dynamics JH, 33] under the 
action of time-dependent magnetic fields and currents. 

Taking H = d Zo E = in Eq. © one recovers the DW dy- 
namics equations for a translationally invariant nanowire with 
no magnetic field applied 12511 . In this case the DW moves 
with a constant velocity (3j/a for small currents j, whereas 
above the critical current j c — na/ \ a — (3 1 in addition to mov- 
ing along the wire its angle rotates around the wire axis. 

Memory device. As an example of the use of equations © 
and © we consider a magnetic memory device based on a flat 
hourglass-shaped nanostrip (shown in Fig. [TJ with the con- 
stant thickness h. We propose a nonvolatile device, which 
employs the magnetization direction within the DW as the in- 
formation storage H44I1 . Without current, the transverse DW 
stays at the place where the nanostrip's cross-section is the 
narrowest. When a particular current pulse is applied, the DW 
magnetization angle <j> flips from to ir. At the intermediate 
step of this flipping process, the DW also deviates from the 
narrowest position of the nanostrip and at the end it comes 
back but with the opposite magnetization direction (for com- 
putational details see Fig. [3] below). The same current pulse 
at a later time can move it back to the original configuration 

ii3. 

The time it takes to switch the magnetization depends on 
the current pulse shape. During this process the main energy 
loss in a realistic wire is the Ohmic loss. How much energy is 
needed for a single switching also depends on the parameters 
of the current pulse. 

We, thus, aim to solve the following problems: 1) What 
is the optimal (requiring the least amount of energy) current 
pulse shape for a given switching time? 2) How the minimal 
required energy per flip depends on the switching time? The 
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answers to these questions are given in Figs. [2] and [3] Here we 
sketch the calculation. 

For a smooth nanostrip, we can approximate its width by a 
parabolic shape w(z) = wo + z 2 /R and its thickness h <C 
Wo, see Fig. [T] Then the correction to the DW energy due 
to the shift from the center of the strip becomes Eq(zq) = 

Jo J 



10000 



Mw AR *o> where J is the exchange constant, A is the DW 
width, and R is the curvature radius of the nanostrip. 

Taking H = (no magnetic field applied) and a zz = 1/ A 
04611 . and rescaling variables as t — > nt/A, zq — > Zo/A, and 
j — > j/n to make all of them dimensionless, Eqs. (0 and © 
become 

j — z + sin 20, (7) 
n(t) = - (a - fi)z + p sin 20 - az = 0. (8) 

Here we have used the fact that in dimensionless variables the 
zq -dependence of the DW energy is Eq(zq) — az 2 with 



7o J A 
MwqkR 




(9) 



FIG. 2. (color online) Dependence of the switching energy in units of 
Jo j 2 M ^ on me switching time T for a = 0.01 (red dashed line), 
a — 0.03 (black solid line) and a = 0.05 (blue dot-dashed line). The 
inset shows the same dependence for smaller range of the switching 
energies. One can clearly see the deviations of the simulated data 
from the analytical solution for the small T limit. 



We can estimate a for materials such as Permalloy using the 
gyromagnetic ratio 70 = 1.76 x 10 11 s _1 T _1 , the exchange 
constant J = 1.3 x 10 11 J/m, the saturation magnetization 
M = 8x 10 5 A/m, and k — j c (a—(3)/a where (a— (3)/ a w 1 
and the critical current measured in units of velocity j c ss 100 
m/s. It gives a m A/(wqR) x 10 nm. Further taking for 
realistic nanowires R = 100 nm, wo = 100 nm, and A ss 10 
nm, we find that a w 0.01 which makes it of the order of a 
and (3. It also justifies neglecting the terms proportional to aa 
and (3a in Eqs. (0 and (©. 

Our goal is to minimize the Ohmic losses IT26I l47Tl per one 
switch of the memory bit. It corresponds to the flipping of 
the DW angle <f> between two stable values to tt, which are 
defined by the minima of transverse anisotropy. To achieve 
this, one has to minimize j 2 during the switching time T while 
keeping the constraint (0. 

The fact that the DW is at rest at the angle = before the 
current pulse and at rest again at the angle <f> — tt immediately 
after the current pulse is taken into account by the boundary 
conditions 0(0) = 0, z (0) = 0, i (+0) = 0, and <j>(T) = tt, 
z (T)=0, io(T-O) =0. 

In order to find minimum of the power J Q T j 2 dt with con- 
straint (0, which must hold at all times, we use the Lagrange 
multiplier J(i) and minimize the functional 



j(t)n(t)]dt 



(in + sin 20) - j(t)n(t) 



dt 



with respect to three functions Zo(t), 4>(t), and J(t) with 
J(0) = 0. Then, in addition to Eq. (©, we obtain two equa- 
tions 

J - 2/3 j cos 20 + 4(i + sin 20) cos 20 = 0, (10) 
i=io + fj-^j, (ID 



where to write the last equation we used Eq. (0, and jo is an 
integration constant. It can be seen that jo = j(0) if J(0) = 0. 

First, we consider the case when the current is absent, j = 
zq + sin 8 = 0. Then J = and equation ([8]) gives 



2a0 cos 20 + a sin 20 = 0. 



(12) 



For small angles 0, this equation describes a damped har- 
monic oscillator. For a <C \[o the motion is underdamped 
and the DW performs oscillations with the frequency ljq = 
\/2o. For a 3> \fo the motion is overdamped. 

If one kicks the DW with a very narrow current pulse j (t) = 
A8(t), integrating io + sin 20 = A6(t) from — e to e with the 
initial condition zq(— e) = and taking the limit e — > 0, we 
find z + = lim z (e) = A. The same procedure for equation 

© gives + = (a — (3)z + — (a — ft) A. After this pulse, the 
motion will be described by equation ( fT2b . Depending on the 
characteristic time scales in this equation, l/y/a or 1/a, the 
motion will be underdamped or overdamped. 

Let us consider the limiting case a = (3 = of the system 
without dissipation. Using Eqs. ©, (©, and ( fTTT l we find 



J = - - + sin 20 - j 
a \ a 



(13) 
(14) 



Substituting these equations into Eq. ( fTOb we also obtain 



d 2 \ ( d 2 
2acos20+— J (asin20+— 0) = 0. (15) 



There is a trivial solution of this equation, 0(i) = 0, but it does 
not satisfy the boundary conditions: the angle does not flip in 
this case. However, in the limit when the switching process is 
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fast, or alternatively a is small, one can find an approximate 
solution. If the switching time T satisfies Ty/a <C 1, all the 
terms except for the fourth derivative can be neglected and we 
obtain d A 4>/dt A = 0. Using the boundary conditions for both 
4> and zo one then finds 

<t>{t) = 3vr(i/T) 2 - 2n(t/T) 3 (16) 

which gives 

j(t) = ±0 + sm20«^=^(l-^y (17) 

The switching energy measured in dimensionless units 
Jq j 2 (t)dt is shown in Fig. [2] For small T, according to 
Eq. (IPTl . it is given by 127t 2 /(ct 2 T 3 ). The actual units of 
the switching energy can be estimated as pJ^A/(j c h), where 
p ps 1.4 X 10~ 7 Sim is Permalloy resistivity, the critical current 
density J c ps 10 12 A/m 2 , A ps h pa 10 nm; and the units of 
time are A/j c ps 10~ 10 s. This gives the units of the switch- 
ing energy to be pa 10~ 16 J. For comparison, best MRAM 
devices typically consume 10~ 10 J per switching a bit [48] in 
5 ns. This shows that our proposed memory device is about 
10 5 times more energy efficient for the same switching times. 
Figure [2] also gives the estimate of the best switching time 
achievable for any particular energy supplied per switch. 

The full set of equations with realistic a and f3 parameters 
can be solved numerically. The comparison of the analytical 
solution (fTTT i for small T with the numerics for a range of T 
is shown in Fig. [3] We find that the linear current coincides 
with the simulation result for the short switching times, i.e. 
T < \j yfd. F° r l° n g switching times, the optimal current 
j (t) flips the DW and then lets it relax with some oscillations 
to zo = and <fi = n within the required time T. 

Summary. We have generalized the DW dynamics to the 
case of translationally non-invariant ferromagnetic nanowires. 
The obtained equations of motion make the description of 
the DW propagation closer to experiment by accounting for 
smooth surface roughness and other disorder effects. We have 
also considered an hourglass-shaped nanostrip with a trans- 
verse DW trapped in the middle as a prototype of a magnetic 
memory device. The exact time-dependence of the current 
pulses required to switch the magnetization with the minimal 
Ohmic losses per switching has been obtained. Furthermore, 
we find how the switching time depends on the Ohmic losses 
per switching for the optimal current pulse. Our estimates 
show that this hourglass-shaped nanodevice may be 10 5 times 
more energy efficient for the same switching times as used in 
the best modern memory devices. 

This work was supported by the NSF Grant No. 0757992, 
ONR-N000141 1 10780, and Welch Foundation (A-1678). 
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FIG. 3. (color online) The time dependence of the DW position and 
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